How Do Differences in Electronic Structure Affect the Use of Vanadium Intermediates as Mimics in Nonheme Iron Hydroxylases?

We study active-site models of nonheme iron hydroxylases and their vanadium-based mimics using density functional theory to determine if vanadyl is a faithful structural mimic. We identify crucial structural and energetic differences between ferryl and vanadyl isomers owing to the differences in their ground electronic states, i.e., high spin (HS) for Fe and low spin (LS) for V. For the succinate cofactor bound to the ferryl intermediate, we predict facile interconversion between monodentate and bidentate coordination isomers for ferryl species but difficult rearrangement for vanadyl mimics. We study isomerization of the oxo intermediate between axial and equatorial positions and find the ferryl potential energy surface to be characterized by a large barrier of ca. 10 kcal/mol that is completely absent for the vanadyl mimic. This analysis reveals even starker contrasts between Fe and V in hydroxylases than those observed for this metal substitution in nonheme halogenases. Analysis of the relative bond strengths of coordinating carboxylate ligands for Fe and V reveals that all of the ligands show stronger binding to V than Fe owing to the LS ground state of V in contrast to the HS ground state of Fe, highlighting the limitations of vanadyl mimics of native nonheme iron hydroxylases.

Machine Learning Prediction of the Experimental Transition Temperature of Fe(II) Spin-Crossover Complexes.

Spin-crossover (SCO) complexes are materials that exhibit changes in the spin state in response to external stimuli, with potential applications in molecular electronics. It is challenging to know a priori how to design ligands to achieve the delicate balance of entropic and enthalpic contributions needed to tailor a transition temperature close to room temperature. We leverage the SCO complexes from the previously curated SCO-95 data set [Vennelakanti et al. J. Chem. Phys. 159, 024120 (2023)] to train three machine learning (ML) models for transition temperature (T1/2) prediction using graph-based revised autocorrelations as features. We perform feature selection using random forest-ranked recursive feature addition (RF-RFA) to identify the features essential to model transferability. Of the ML models considered, the full feature set RF and recursive feature addition RF models perform best, achieving moderate correlation to experimental T1/2 values. We then compare ML T1/2 predictions to those from three previously identified best-performing density functional approximations (DFAs) which accurately predict SCO behavior across SCO-95, finding that the ML models predict T1/2 more accurately than the best-performing DFAs. In addition, we study ML model predictions for a set of 18 SCO complexes for which only estimated T1/2 values are available. Upon excluding outliers from this set, the RF-RFA RF model shows a strong correlation to estimated T1/2 values with a Pearson's r of 0.82. In contrast, DFA-predicted T1/2 values have large errors and show no correlation to estimated T1/2 values over the same set of complexes. Overall, our study demonstrates slightly superior performance of ML models in comparison with some of the best-performing DFAs, and we expect ML models to improve further as larger data sets of SCO complexes are curated and become available for model training.

Why Nonheme Iron Halogenases Do Not Fluorinate C-H Bonds: A Computational Investigation.

Selective halogenation is necessary for a range of fine chemical applications, including the development of therapeutic drugs. While synthetic processes to achieve C-H halogenation require harsh conditions, enzymes such as nonheme iron halogenases carry out some types of C-H halogenation, i.e., chlorination or bromination, with ease, while others, i.e., fluorination, have never been observed in natural or engineered nonheme iron enzymes. Using density functional theory and correlated wave function theory, we investigate the differences in structural and energetic preferences of the smaller fluoride and the larger chloride or bromide intermediates throughout the catalytic cycle. Although we find that the energetics of rate-limiting hydrogen atom transfer are not strongly impacted by fluoride substitution, the higher barriers observed during the radical rebound reaction for fluoride relative to chloride and bromide contribute to the difficulty of C-H fluorination. We also investigate the possibility of isomerization playing a role in differences in reaction selectivity, and our calculations reveal crucial differences in terms of isomer energetics of the key ferryl intermediate between fluoride and chloride/bromide intermediates. While formation of monodentate isomers believed to be involved in selective catalysis is shown for chloride and bromide intermediates, we find that formation of the fluoride monodentate intermediate is not possible in our calculations, which lack additional stabilizing interactions with the greater protein environment. Furthermore, the shorter Fe-F bonds are found to increase isomerization reaction barriers, suggesting that incorporation of residues that form a halogen bond with F and elongate Fe-F bonds could make selective C-H fluorination possible in nonheme iron halogenases. Our work highlights the differences between the fluoride and chloride/bromide intermediates and suggests potential steps toward engineering nonheme iron halogenases to enable selective C-H fluorination.

Assessing the performance of approximate density functional theory on 95 experimentally characterized Fe(II) spin crossover complexes.

Spin crossover (SCO) complexes, which exhibit changes in spin state in response to external stimuli, have applications in molecular electronics and are challenging materials for computational design. We curate a dataset of 95 Fe(II) SCO complexes (SCO-95) from the Cambridge Structural Database that have available low- and high-temperature crystal structures and, in most cases, confirmed experimental spin transition temperatures (T1/2). We study these complexes using density functional theory (DFT) with 30 functionals spanning across multiple rungs of "Jacob's ladder" to understand the effect of exchange-correlation functional on electronic and Gibbs free energies associated with spin crossover. We specifically assess the effect of varying the Hartree-Fock exchange fraction (aHF) in structures and properties within the B3LYP family of functionals. We identify three best-performing functionals, a modified version of B3LYP (aHF = 0.10), M06-L, and TPSSh, that accurately predict SCO behavior for the majority of the complexes. While M06-L performs well, MN15-L, a more recently developed Minnesota functional, fails to predict SCO behavior for all complexes, which could be the result of differences in datasets used for parametrization of M06-L and MN15-L and also the increased number of parameters for MN15-L. Contrary to observations from prior studies, double-hybrids with higher aHF values are found to strongly stabilize high-spin states and therefore exhibit poor performance in predicting SCO behavior. Computationally predicted T1/2 values are consistent among the three functionals but show limited correlation to experimentally reported T1/2 values. These failures are attributed to the lack of crystal packing effects and counter-anions in the DFT calculations that would be needed to account for phenomena such as hysteresis and two-step SCO behavior. The SCO-95 set thus presents opportunities for method development, both in terms of increasing model complexity and method fidelity.

Non-Native Anionic Ligand Binding and Reactivity in Engineered Variants of the Fe(II)- and α-Ketoglutarate-Dependent Oxygenase, SadA.

Mononuclear non-heme Fe(II)- and α-ketoglutarate-dependent oxygenases (FeDOs) catalyze a site-selective C-H hydroxylation. Variants of these enzymes in which a conserved Asp/Glu residue in the Fe(II)-binding facial triad is replaced by Ala/Gly can, in some cases, bind various anionic ligands and catalyze non-native chlorination and bromination reactions. In this study, we explore the binding of different anions to an FeDO facial triad variant, SadX, and the effects of that binding on HO• vs X• rebound. We establish not only that chloride and bromide enable non-native halogenation reactions but also that all anions investigated, including azide, cyanate, formate, and fluoride, significantly accelerate and influence the site selectivity of SadX hydroxylation catalysis. Azide and cyanate also lead to the formation of products resulting from N3•, NCO•, and OCN• rebound. While fluoride rebound is not observed, the rate acceleration provided by this ligand leads us to calculate barriers for HO• and F• rebound from a putative Fe(III)(OH)(F) intermediate. These calculations suggest that the lack of fluorination is due to the relative barriers of the HO• and F• rebound transition states rather than an inaccessible barrier for F• rebound. Together, these results improve our understanding of the FeDO facial triad variant tolerance of different anionic ligands, their ability to promote rebound involving these ligands, and inherent rebound preferences relative to HO• that will aid efforts to develop non-native catalysis using these enzymes.

Probing the Mechanism of Isonitrile Formation by a Non-Heme Iron(II)-Dependent Oxidase/Decarboxylase.

The isonitrile moiety is an electron-rich functionality that decorates various bioactive natural products isolated from diverse kingdoms of life. Isonitrile biosynthesis was restricted for over a decade to isonitrile synthases, a family of enzymes catalyzing a condensation reaction between l-Trp/l-Tyr and ribulose-5-phosphate. The discovery of ScoE, a non-heme iron(II) and α-ketoglutarate-dependent dioxygenase, demonstrated an alternative pathway employed by nature for isonitrile installation. Biochemical, crystallographic, and computational investigations of ScoE have previously been reported, yet the isonitrile formation mechanism remains obscure. In the present work, we employed in vitro biochemistry, chemical synthesis, spectroscopy techniques, and computational simulations that enabled us to propose a plausible molecular mechanism for isonitrile formation. Our findings demonstrate that the ScoE reaction initiates with C5 hydroxylation of (R)-3-((carboxymethyl)amino)butanoic acid to generate 1, which undergoes dehydration, presumably mediated by Tyr96 to synthesize 2 in a trans configuration. (R)-3-isocyanobutanoic acid is finally generated through radical-based decarboxylation of 2, instead of the common hydroxylation pathway employed by this enzyme superfamily.

Harder, better, faster, stronger: Large-scale QM and QM/MM for predictive modeling in enzymes and proteins.

Computational prediction of enzyme mechanism and protein function requires accurate physics-based models and suitable sampling. We discuss recent advances in large-scale quantum mechanical (QM) modeling of biochemical systems that have reduced the cost of high-accuracy models. Tradeoffs between sampling and accuracy have motivated modeling with molecular mechanics (MM) in a multiscale QM/MM or iterative approach. Limitations to both conventional density-functional theory and classical MM force fields remain for describing noncovalent interactions in comparison to experiment or wavefunction theory. Because predictions of enzyme action (i.e. electrostatics), free energy barriers, and mechanisms are sensitive to the protocol and embedding method in QM/MM, convergence tests and systematic methods for quantifying QM-level interactions are a needed, active area of development.

Molecular basis of C-S bond cleavage in the glycyl radical enzyme isethionate sulfite-lyase.

Desulfonation of isethionate by the bacterial glycyl radical enzyme (GRE) isethionate sulfite-lyase (IslA) generates sulfite, a substrate for respiration that in turn produces the disease-associated metabolite hydrogen sulfide. Here, we present a 2.7 Å resolution X-ray structure of wild-type IslA from Bilophila wadsworthia with isethionate bound. In comparison with other GREs, alternate positioning of the active site ß strands allows for distinct residue positions to contribute to substrate binding. These structural differences, combined with sequence variations, create a highly tailored active site for the binding of the negatively charged isethionate substrate. Through the kinetic analysis of 14 IslA variants and computational analyses, we probe the mechanism by which radical chemistry is used for C-S bond cleavage. This work further elucidates the structural basis of chemistry within the GRE superfamily and will inform structure-based inhibitor design of IsIA and thus of microbial hydrogen sulfide production.

When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins.

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O-H⋯O HBs than N-H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction.